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We extend the real-space renormalization group (RG) approach to the study of the energy level 
statistics at the integer quantum Hall (QH) transition. Previously it was demonstrated that the RG 
approach reproduces the critical distribution of the po«;er transmission coefficients, i.e., two-terminal 
conductances, Pc{G), with very high accuracy. The RG flow of P{G) at energies away from the 
transition yielded the value of the critical exponent, that agreed with most accurate large-size 
lattice simulations. To obtain the information about the level statistics from the RG approach, 
we analyze the evolution of the distribution of phases of the amplitude transmission coefhcient 
upon a step of the RG transformation. From the fixed point of this transformation we extract the 
critical level spacing distribution (LSD). This distribution is close, but distinctively different from 
the earlier large-scale simulations. We find that away from the transition the LSD crosses over 
towards the Poisson distribution. Studying the change of the LSD around the QH transition, we 
check that it indeed obeys scaling behavior. This enables us to use the alternative approach to 
extracting the critical exponent, based on the LSD, and to find v = 2.37 ± 0.02 very close to the 
value established in the literature. This provides additional evidence for the surprising fact that a 
small RG unit, containing only five nodes, accurately captures most of the correlations responsible 
for the localization-delocalization transition. 

PACS numbers: 73.43.-f, 73.43.Nq, 64.60.Ak 



I. INTRODUCTION 

It has been realized long ago that, alongside with 
the change in the behavior of the eigcnfunctions, a 
localization-delocalization transition manifests itself in 
the statistics of the energy levels. In particular, as the 
energy is swept across the mobility edge, the shape of 
the level spacing distribution (LSD) crosses over from 
the Wigner-Dyson distribution, corresponding to the ap- 
propriate universality class, to the Poisson distribution. 

Moreover, finite-size corrections to the critical LSD ex- 
actly at the mobility edge allow tp, determine the value 
of the correlation length exponent p thus avoiding an ac- 
tual analysis of the spatial extent of the wave functions. 
For this reason, the energy level statisticSj^cnnstitutes 
an alternative to the MacKinnon-KrameiBdu^El and to 
the transmission-matrixBQ approaches to the numerical 
study of localization. 

Another reason why a large number of 
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ried out during the past decade is the controversy 



,±ail of the critical 
JEI that this tail is 



that existed over the large-spa 
LSD. Conclusive demonstratior£3ll^ 
Poissonian, i.e., that there is no repulsion between the. 
levels with spacings much larger than the mean value,Ej 
rather than super-Poissonian,Ej implying that repulsion 
is partially preseriied- required a very high accuracy 
of the simulations Ej'EI The bulk of numerical work on 
the level statistics at the transition was carried out 



for three-dimensional systemsB§00lli'0'0 for which 
there exists a mobility edge separating localized and 
extended states. In two dimensions all the states are 
localized in the absence of a magnetic field. In the 
presence of a magnetic field, localization-delocalization 
transitions in two dimensions (quantum Hall transitions) 
arc infinitely sharp. Still the reasoning of Ref. |l| applies. 
Numerical studies have established a Poissonian tail 
of the LSD .113 It was also demonstratectj that the 
procedure of extracting the localization length exponent 
from the finite-size corrections yields a value close to 
v = 2.35 Jound .from large-size simulations of the wave 
functions EflMBEZi 

Recently, a semianalytical description of the inte- 
ger quantum Hall transition, based on the extepsion 
of the scaling ideas for the classical percolations to 
the Chalkep-Coddington (CC) modeLof the quantum 
percolationjH3 has been developed Hd'Eil The key idea of 
this description, a real-space-renormalization group ap- 
proach (RG), is the following. Each RG step corresponds 
to a doubling of the system size. The RG transformation 
relates the conductance distribution of the sample at the 
next step to the conductance distribution at the previous 
step. The fixed point of this transformation, yields the 
distribution of the conductance, Pc{G) of a macroscopic 
sample at the quantum Hall transition. This universal 
distribution describes the mesoscopic properties of a fully 
coherent quantuna-Hall sample. Analogously to the clas- 
sical percolation,L3 the correlation length exponent, v, 
was extracted from the RG procedureEHl using the fact 
that a slight shift of the initial distribution with respect 
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to the fixed-point, Pc{G), drives the system to the insu- 
lator upon renormahzation. Then the rate of the shift of 
the distribution maximum determines the value of z/. Re- 
markably, both Pc{G) an d . f he c ritical exponent obtained 
within the RG approachUaEla agree very wellpariiJi^he 
"exact" results of the large-scale simulations .l§EZIHc3'LZI 
The goal of the present paper is two-fold. Firstly, we 
extend the RG approach to the level statistics at the 
transition in order to subject its validity to yet another 
test. Secondly, we apply the method analogous to the 
finite-size-corrections analysis to extract u from the LSD 
obtained within the RG approach. This method yields 
V = 2.37 ± 0.02, which is even closer to the most pre- 
cise large-scale simulations result i> = 2.35 ± 0.03E3 than 
the value v :p-|2.39 ± 0.01 inferred from the conductance 
distribution.Ea The latter result is by jap means trivial. 
Indeed, the original RG transformations related the con- 
ductances, i.e., the absolute values of the transmission co- 
efficients of the original and the doubled samples, while 
the phases of the transmission coefficients were assumed 
random and uncorrelated. In contrast, the level statis- 
tics at the transition corresponds to the fixed point in 
the distribution of these phases. Therefore, the success 
of the RG approach for conductances does not guarantee 
that it will be equally accurate quantitatively for the level 
statistics. 

Within both RG transformations, for the magnitudes 
and for the phases of the transmission coefficients, an 
initial deviation from the critical distribution drives the 
system towards an insulator with zero transmission and 
Poissonian LSD. Thus the procedures of the extraction 
of v from both transformations are technically differ- 
ent, but conceptually similar. In fact, the shape of the 
critical LSD, obtained from the RG approach, shows 
systeniatic,. ck'V ii ^yo ff^ r^*^"^ large-scale simulation 
resultsOMBEJElMEaH which yield the body of LSD 
very close tp-thc Gaussian unitary random matrix ensem- 
ble (GUE).N However, the RG flow of the LSD towards 
the insulator appears to be robust. 

The paper is organized as follows iFiKSt, in Sec. || we 
review the real-space RG approachEjEiftJ and adjust it 
to the computation of the energy levels and the LSD. 
In Sec. [II we present our numerical results for the LSD. 
The finite-size scaling (ESS) analysis of the obtained LSD 
at the QH transition is reported in Sec. IV . Concluding 
remarks are presented in Sec. 



based on the RG unit shown in Fig. |^. The unit is a frag- 
ment of the CC network consisting of five nodes. Each 
node, i, is characterized by the transmission coefficient 
ti, which is an amplitude to deflect an incoming electron 




FIG. 1: Chalker-Coddington network on a square lattice con- 
sisting of nodes (circles) and links (arrows). The RG unit 
used for Eq. (^ combines five nodes (full circles) by neglecting 
some connectivity (dashed circles). $i, . . . , $4 are the phases 
acquired by an electron along the loops as indicated by the 
arrows, ^l/i, . . . , »I'4 represent wave function amplitudes, and 
the thin dashed lines illustrate the boundary conditions used 
for the computation of level statistics. 



II. MODEL AND RG METHOD FOR THE LSD 

A. RG approach to the conductance distribution 

A detailed description of the RG approach to the con- 
ductance distribution can be found in Refs. |30|,|3l|,P4 It is 



along the link to the left. Analogously, the reflection co- 
efficient = (1 — tfy/^ is the amplitude to deflect the 
incoming electron to the right. Doubling of the sample 
size corresponds to the replacement of the RG unit by a 
single node. The RG transformation expresses the trans- 
mission coefflcient of this effective node, t' , through the 
transmission coefficients of the five constituting nodea23 
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t' = 



tMr2r3rie"^^ - 1) + t2t4e'(*='+**Hnr3r5e-'*i - 1) + ^3(^2^56'*^ + ^1*46'**) 



I — 



(1) 



Here $j are the phases accuniulated along the closed 
loops (see Fig. |l|). Within the RG approach to the 
conductance distribution, information about electron en- 
ergy is incorporated only into the values of ti. The en- 
ergy dependence of phases, is irrelevant; they are 
assumed completely random. Due to this randomness, 
the transmission coefficients, ti, for a given energy, are 
also randomly distributed with a distribution function 
P{t). Then the transformation (|l|) allows, upon averag- 
ing over <i>j, to generate the next-step distribution P{t'). 
Therefore, within the RG scheme, a delocalized state cor- 
responds to the fixed point, Pc{t), of the RG transforma- 
tion. Due to the symmetry of the RG unit, it is obvious 
that the critical distribution, Pc(i^), of the power trans- 
mission coefficient, = G, which has the meaning of 
the two-terminal conductance, is symmetric with respect 
to = i. In other words, the RG transformation re- 
spects the duality between transmission and reflection. 
The critical distribution Pc(G) found in Refs. |3^ and ^ 
agrees very well with the results of direct large-scale sim- 
ulations. 



B. RG approach to the LSD 

Universal features of the energy level statistics in a 
macroscopic fully coherent sample at the quantum Hall 
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transition complement the universality in the conduc- 
tance distribution. The prime characteristics of the level 
statistics is the LSD - the distribution of the spacings be- 
tween neighboring energy levels. In order to adjust the 
RG approach to the calculation of the LSD, it is neces- 
sary to "close" the sample at each RG step in order to 
discretize the energy levels. One of the possible variants 
of such a closing is shown in Fig. |l| with dashed lines. 

For a given closed RG unit with a fixed set of t^-values 
at the nodes, the positions of the energy levels are de- 
termined by the energy dependences, ^j{E), of the four 
phases along the loops. These phases change by ~ tt 
within a very narrow energy interval, inversely propor- 
tional to the sample size. Within this interval the change 
of the transmission coefficients is negligibly small. A 
closed RG unit in Fig. |^ contains 10 links, and, thus, 
it is described by 10 amplitudes. These amplitudes are 
related by 10 equations (2 at each node). Each link is 
characterized by an individual phase. On the other hand, 
it is obvious that the energy levels are determined only 
by the phases along the loops. One possible way to derive 
the system, in which individual phases combine into $j 
is to exclude from the original system of 10 equations all 
amplitudes except the "boundary" amplitudes (see 
Fig. 1). This procedure is similar to the derivation of 
Eq. (1). The system of equations for the remaining four 
amplitudes takes the form 
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where the parameter lo should be set zero. Then the 
energy levels, Ek, of the closed RG unit are the energies 
for which, with phases $j (E) — [E^. ) , one of the four 
eigenvalues of the matrix in the left-hand side of Eq. (||) 
is equal to one. If we keep uj in the right-hand side of 
Eq. (^, then the above condition can be reformulated 
as <jj{Ek) ~ 0. Thus, the calculation of the energy levels 
reduces to a diagonalization of the 4x4 matrix. 

The crucial step now is the choice of the dependence 
^j{E). If each loop in Fig. |l| is viewed as a closed 
equipotentiaJ as it is the case for the first step of the RG 
procedure,E3 then ^j{E) is a true magnetic phase, which 
changes linearly with energy with a slope governed by 



the actual potential profile, which, in turn, determines 
the drift velocity. Thus we have 



<S>j{E) = <S>oj+2tt 



E 



(3) 



where a random part, ^o.j, is uniformly distributed 
within [0,27r], and ^ir/sj is a random slope. Strictly 
speaking, the dependence (^) applies only for the first 
RG step. At each following step, n > 1, ^j{E) is a com- 
plicated function of E which carries information about 
all energy scales at previous steps. However, in the spirit 
of the RG approach, we assume that ^j{E) can still be 
linearized within a relevant energy interval. The conven- 
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tional RG approach suggests that different scales in the 
real space can be decoupled. Linearization of Eq. (^) im- 
plies a similar decoupling in the energy space. In the case 
of phases, a "justification" of such a decoupling is that at 
each following RG step, the relevant energy scale, that is 
the mean level spacing, reduces by a factor of 4. 

With ^j{E) given by Eq. (||) and fixed values of ij, 
the statistics of energy levels determined by the matrix 
equation ^ is obtained by averaging over the random 
initial phases <i>oj - In particular, each realization of $oj 
yields 3 level spacings which are then used to construct 
a smooth LSD. We now outline the RG procedure for 
the LSD. The slopes Sj in Eq. (||) determine the level 
spacings at the first step. They are randomly distributed 
with a distribution function Pq{s). Diagonalization of the 
matrix in Eq. (j^) with subsequent averaging over realiza- 
tions yields the LSD, i-i(s), at the second step. Then the 
key element of the RG procedure, as applied to the level 
statistics, is using Pi{s) as a distribution of slopes in Eq. 
(^. This leads to the next-step LSD and so on. 

It is instructive to compare our procedure of calcu- 
lating the energy levels with an approach adopted in 
large-scale simulations within the CC modeLO'Ej This, 
approach is based on the unitary network operator Ux3 
For a single RG unit this operator acts analogously to 
the matrix in the left-hand side of Eq. (||). However, 
within the approach of Refs. ^and the energy depen- 
dence of phases in the elements of the matrix was ne- 
glected (only the random contributions, $oji were kept). 
Then, instead of the energy levels, Ek, diagonalization 
of the matrix (^) yielded a set of eigenvalues, exp(iwfe). 
The numbers ujk were named quasienergies, and it is the 
statistics of these quasienergies that was studied in Ref. 
|39| . Comparison of the two procedures for a single RG 
unit is illustrated in Fig. ||. Fig. || shows the dependence 
of the 4 quasienergies ujk on the energy E calculated for 
two single sample RG units, with ti chosen from the crit- 
ical distribution Pc{t). The energy dependence of the 
phases was chosen from LSD of the GUE according 
to Eq. (H). It is seen that the dependences uj{E) range 
from remarkably linear and almost parallel (Fig. ^ja) to 
strongly nonlinear (Fig. ^). 
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FIG. 2: Energy dependence of the quasieigenenergies uj for 
two sample configurations. Instead of using the quasispec- 
trum obtained from uui{E = 0) (O) we calculate the real 
eigenenergies according to (jj{Ek) — (□). Diflterent line styles 
distinguish different uii{E). We emphasize that the observed 
behavior varies from sample to sample between remarkably 
linear (a) and strongly nonlinear (b). 



III. NUMERICAL RESULTS 
A. The LSD at the QH transition 

As a first step of the RG procedure for the calculation 
of the critical LSD we chose for Pa{s) the distribution cor- 
responding to the GU.F .random matrix ensemble, since 
previous simulationalBo indicated that the LSD at the 
transition is close to GUE. According to Po{s), we pick 
Sj and set <i>j, j = 1, . . . , 4 as in Eq. (^). For the trans- 
mission coefficients_ii , i — 1, . . . , 5 we usej_the fixed-point 
distribution P^ {t) ,cd obtained previously.LJ 

From the solutions of Eq. (^) corresponding to 
ujj{Ek) = the new LSD Pi{s') is constructed using the 



"unfolded" energy level spacings sj„ — (£^„i+i — E,n)/A, 
where m = 1,2,3, Ek+i > Ek and the -mean spacing 
A = {E4- El)/?,. Due to the "unfolding"El with A, the 
average spacing is set to one for each sample and in each 
RG-iteration step we superimpose spacing data of 2 x 10^ 
RG units. The resulting LSD is discretized in bins with 
largest width 0.01. In the following iteration step we re- 
peat the procedure using Pi as initial distribution. We 
assume that the iteration process has converged when the 
mean-square deviation of distribution Pn{s) deviates by 
less than 10"'* from its predecessor P„_i(s). The RG it- 
eration process converges rather quickly after only 2 — 3 
RG steps. The resulting LSD, Pc(s), is shown in Fig. 
^ together with an LSD for the unitary random matrix 




FIG. 3: Critical distributions Pc(s) obtained from the spec- 
trum of uJi{E — 0) and from the RG approach using the real 
eigenenergies Ek in comparison to the LSD for GUE. As in all 
other graphs P{s) is shown in units of the mean level spacing 
A. 



ensemble. 

Although Pc{s) exhibits the expected features, namely, 
level repulsion for small s and a long tail at large s, the 
overall shape of Pc{s) differs noticeabljj^pin GUE. In the 
previous large-size lattice simulationsllSEj the obtained 
critical LSD was much closer to GUE than Pc{s) in Fig. 
^. This fact, however, does not reflect on the accuracy 
of the RG approach. Indeed, as it was demonstrated 
recently, the critical LSD - although being system size 
independent — p-nevertheless depends on the geometry 



of the samples 2 
ary conditions .cf 



-ajid on the specific choice of bound- 
O Sensitivity to the boundary condi- 
tions does not affect the asymptotics of the critical dis- 
tribution, but rather manifests itself in the shape of the 
"body" of the LSD. Recall now that the boundary con- 
ditions which we have imposed to calculate the energy 
levels (dashed lines in Fig. |^) are non-periodic. 

There is another possibility to assess the critical LSD, 
namely by iterating the distribution of quasienergies. In 
Fig. we show the result of this procedure. It appears 
that the resulting distribution is almost identical to Pc{s)- 
This observation is highly non-trivial, since, as follows 
from Fig. |^, there is no simple relation between the ener- 
gies and quasienergies. Moreover, if instead of the linear 
i?-dependence of $j, we choose another functional form, 
say. 



$j(£') = + 2arcsin 



E 



2p 



(4) 



2p 



< 1, then, the 



where the integer p insures that 

RG procedure would yield an LSD which is markedly dif- 
ferent (within the "body") from Pc{s). This is illustrated 
in Fig. ^. 



FIG. 4: Critical distributions Pc{s) for a linear and an arcsin 
energy dependence of the phases $j. The form of Pc{s) is 
clearly influenced by the actual choice of ^j{E). Hence the 
balk of the distribution is non-universal. The inset illustrates 
examples of the two different functions ^j{E) as in Eqs. (^) 
and (0). 



Both procedures, using quasienergies instead of real 
and linearization of the energy 



energies [as in Ref. |39[ 

dependence of phases [as in Eq. (||)] are not rigorous 
Linearization is dictated by the RG concept. The coinci- 
dence of the results of the two procedures indicates that 
the concept of quasienergies, namely, that they obey the 
same statistics as real energies, is equivalent to the RG. 



B. Small and large s behavior 

As it was mentioned above, the general shape of the 
critical LSD is not universal. However, the small s behav- 
ior of Pc{s) must be the same as for the unitary random 
matrix ensemble, namely Pc{s) oc s^. This is because 
delocalization aJjfil|ie quantum Hall transition implies the 
level repulsijOal 
critical LSDEIB 




simulations of the 
B satisfy this general 
requirement. The same holds also for our result, as can 
be seen in Fig. ^. The given error bars of our numerical 
data are standard deviations computed from a statistical 
average of 100 FP distributions each obtained for differ- 
ent random sets of t^'s and $j's within the RG miit. In 
general, within the RG approach, the s^-asymptotics of 
P{s) is most natural. This is because the levels are found 
from diagonalization of the 4x4 unitary matrix with ab- 
solute values of elements widely distributed between 
and 1. 

The right form of the large-s tail of P(s) is Poissonian, 
Pc(s) C!C exp(— 6s).Ej For the Anderson model in three di- 
mensions, unambiguous confirmation of this prediction 
in numerical simulations became possible onb/ .when very 
high numerical accuracy had been achieved. E3il3 This is 
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FIG. 5: Critical Pc{s) for small s in agreement with the pre- 
dicted s'^ behavior. Due to the log-log plot errors are shown 
in the upper direction only. 

because Pc{s) assumes the Poissonian asymtotics only at 
large enough s ^ 3A. For the quantum Hall transition, 
a linear behavior of InPc(s) with a slope corresponding 
to the value 6 « 4.1 has been found in Ref. |l^ from the 
analysis of the interval 2 < s/A < 4. Our data, as shown 
in Fig. ^, has a high accuracy only for s/A < 2.5. For 
such s, the distribution Pc{s) does not yet reach its large- 
s tail. Thus, the value of parameter b extracted from this 
limited interval is somewhat ambiguous. Namely, we ob- 
tain b = 5.442 for s/A e [1.5,2.0] and b = 6.803 for 
s/A e [2.0,2.5]. 

Summarizing, the accuracy of the RG approach, ap- 
plied to the level statistics, is insufficient to discern the 
only non-trivial feature of the critical LSD, i.e., the uni- 
versal Poissonian asymtotics. However, the scaling anal- 
ysis of LSD clearly reveals the universal features of the 
quantum Hall transition as we demonstrate in the next 
Section. 



IV. SCALING RESULTS FOR THE LSD 
A. Finite-size scaling at the QH transition 

The critical exponent, v, of the quantum Hall transi- 
tion governs the divergence of the correlation length 
as a function of the arbitrary control parameter zq, i.e., 

^ooizo) cc\za - Zcl^", (5) 

where Zc is the critical value. The values of v calculated, 
using diffecent numerical methods, e.g., z/ — 2.35±0.03£3 
2.4 ± 0.2P 2.5 ± 0.5EI agree with each other. The RG 
approach for the conductance distribution also yields a 
rather accurate value v — 2.39 ± 0.01.l3 In Sec. II we 
have introduced a complimentary RG approach to the 
distribution of the energy levels at the transition. It can 



FIG. 6: The large s tail of IsS.s) compared with fits according 
to the predictions of Ref. hi (lines). The interval used for 
fitting is indicated by the bars close to the lower axis. For 
clarity errors are shown in upper direction and for s/A = 
1.5,2.0,2.5,3.0 only. For s/A < 2.4, only every 5th data 
point is drawn by a symbol. 

be expected on general grounds, that the LSD obtained 
from the RG approach obeys scaling at small enough zq — 
Zc- However, it is by now means obvious, whether the 
value of V extracted from different variants of the RG 
approach are consistent. 

In order to extract v from the LSD wc employ the 
one-parameter-scaling analysis. This analysis is based 
on the rescaling of a quantity a{N; {zi}) — depending 
on (external) system parameters {zi} and the system size 
TV — onto a single curve by using a scaling function / 

Since Eq. (^), as indicated by "oo" , holds only in the limit 
of infinite system size, we now use the scaling assumption 
to extrapolate / to TV — > oo from the finite-size results of 
the computations. Once / and ^oo are known, the value 
of ly can be then inferred. _ 

In the original formulation of the RG approachEj it 
was demonstrated that there is a natural parametrization 
of the transmission coefficients t, i.e., t — {e^ + 1)^^/^. 
For such a parametrization, z can be identified with a 
dimensionless electron energy. The quantum Hall transi- 
tion occurs at z = 0, which corresponds to the center of 
the Landau band. The universal conductance distribu- 
tion at the transition, Pc{G), corresponds to the distribu- 
tion Qc{z) = Pc [(e^ + 1)^1] /4cosh^(z/2) of parameter 
z, which is symmetric witb— respect to z = O and has 
a shape close to a gaussian.E3 The RG procedure for the 
conductance distribution converges and yields Qc{z) only 
if the initial distribution is an even function of z. This 
suggests to choose as a control parameter in Eq. (||), zq, 
the position of the maximum of the function Q(z). Then 
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the meaning of zq is the electron energy measured from 
the center of Landau band. The fact that the quan- 
tum Hall transition is infinitely sharp implies that for 
any zq 7^ 0, the RG procedure drives the initial distribu- 
tion Q{z — Zq) towards an insulator, either with complete 
transmission of the network nodes (for zq > 0) or with 
complete reflection of the nodes (for zq < 0). 



B. Scaling for ap and ai 

In principle, we are free to choose for the finite-size 
analysis any characteristic quantity a(iV; zq) constructed 
from the LSD which has a systematic dependence on sys- 
tem size N for zq 7^ while being constant at the tran- 
sition Zft-.p.£L|_Bejause of the large number of possible 
choicealla'EjEjcZica we restrict ourselves to two quantities 
which are obtained by integration of the LSD and have 
already been successfully used in Rcfs. || and namely 



ap 



P{s)ds 



and second 



1 f''° 
ai = — / I{s)ds, 
So Jo 



(7) 



(8) 



with I{s) — P{s')ds' . The integration limit is chosen 
as son= 1-4 which approximates the common crossing 
points of all LSD curves as can be seen in Fig. 0. Thus 
P(so) is independent of the distance |zo — to the crit- 
ical point and the system size N. We note that N is 
directly related to the RG step n by = 2". The dou- 
ble integration in ai is numerically advantageous since 
fluctuations in P{s) are smoothed. We now apply the 
finite-size-scaling approach from Eq. (p|) 



ai^p{N,zo) = f 



N 



^ooizo) 



(9) 



Since ai^p{N, zq) is analytical for finite N, one can ex- 
pand the scaling function / at the critical point. The 
first order approximation yields 



a{N, Zq) « aiN, z^) + a\za - z^lN^/" 



(10) 



where a is a dimensionless coefficient. For our calculation 
we use ariigher order expansion proposed by Slevin and 
Ohtsuki.EJ In Ref. ^ the function / is expanded twice, 
first, in terms of the Chebyshev polynomials of order 
and, second, in Taylor series with terms |zo — Zc| in the 
power Oz. This procedure allows to describe the devi- 
ations from linearity in j^g — -^cl at the transition. In 
addition, in Ref. ^ the contributions from an irrelevant 
scaling variable which leads to a shift of the transition 
for small system sizes was taken into account. In our 
case, in contrast to the Anderson model of localization, 
the transition point Zq = is known. Therefore, we can 
neglect the influence of irrelevant variables. In order to 
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FIG. 7; RG of the LSD used for the computation of The 
dotted lines corresponds to the first 9 RG iterations with 
an initial distribution Pq shifted to complete transmission 
{zo — 0.1) while the long-dashed lines represent results for 
a shift toward complete reflection {zq = —0.1). Within the 
RG procedure the LSD moves away from the FP as indicated 
by the arrows. At s/A « 1.4 the curves cross at the same 
point - a feature we exploit when deriving a scaling quantity 
from the LSD. 



obtain the functional form of /, the fitting parameters, 
including v, are evaluated by a nonlinear least-square 
(x^) minimization. In Fig. |^ we show the resulting fit for 
ap and ai at the transition. 

The fits are chosen in such a way that the total num- 
ber of parameters is kept at a minimai-value, while the fit 
agrees well with the numerical data.Eil The correspond- 
ing scaling curves are displayed in Fig. ^. In the plots 
the two branches corresponding to complete reflection 
(zq < 0) and complete transmission (zq > 0) can be 
clearly distinguished. In order to estimate the error of 
fitting procedure we compare the results for ly obtained 
by different orders and Oz of the expansion, system 
sizes A, and regions around the transition. A part of our 
over 100 fit results together with the standard deviation 
of the fit are given in Table |. The value of v is calculated 
as the average of all individual fits where the resulting 
error of f was smaller than 0.02. The error is then de- 
termined as the standard deviation of the contributing 
values. By this method we assure that our result is not 
influenced by local minima of the nonlinear fit. So we 
consider 1/ — 2.37 ± 0.02 as a reliable value for the expo- 
nent of the localization length at the QH transition ob- 
tained from the RG approach to LSD. This is in excellent 
agreement with 1/ = 2.35 ± 0.03 (ReL ||), 2.4 ± 0.2 (ReL 
13), 2.5 ± 0.5 (ReL and 2.39 ± 0.01 (ReL |3|) calcu- 
lated previously. In addition to ap and ai, we tested also 
a parameter-free scaling quantity s^P(s)ds,c3 where 
the whole distribution P{s) is taken into account. Here, 
due to the influence of the large s-tail a less reliable value 
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FIG. 8: Behavior of ai and ap at the QH transition as re- 
sults of the RG of the LSD. Data are shown for RG iter- 
ations n = 1,...,9 corresponding to effective system sizes 
N = 2" = 2,...,512. Full lines indicate the functional de- 
pendence according to FSS using the minimization with 
= 2 and = 3. 



FIG. 9: Finite size scaling curves resulting from the fit of 
our data shown in Fig. H. Different symbols correspond to 
different effective system sizes N = 2". The data points col- 
lapse onto a single curve indicating the validity of the scaling 
approach. 



p = 2.33 ± 0.05 was obtained. 



C. Test of consistency 

Finally we address the question, how the actual form 
of the distribution Q{z) affects the results for LSD and 
the scaling analysis. Recall that in the above calcula- 
tions we have used at each step of the RG procedure the 
distribution Q{z) derived from the critical conductance 
distribution, Pc{G). The function Pc{G) is shown in Fig. 

(inset) with a full line. In order to understand the 
importance of the fact that Pc{G) is almost flat, we have 
repeated our calculations choosing for P{G) a relatively 
narrow gaussian distribution P{G) = PcmEiG) at each 
RG step. This distribution is shown with a dashed line in 
Fig. ^(inset). The obtained LSD is presented in Fig. 
Obviously, it agrees much worse with GUE, which can be 
considered as a reference point, than the LSD computed 
using the true Pc{G). Our data for ai calculated for 



P{G) = PGauB(G) is plotted in Fig. |ri|. The curves for 
small system sizes N exhibit strong deviations, i.e., there 
is initially no common crossing point, while for large N 



a behavior similar to Fig 
N data are neglected in t 



is observed. Therefore, small 
ihe scaling analysis. The fits 
for ai and ap are carried out using zq G [—0.05, 0.05] and 
N ^ 16- 512. They yield the values i/j = 2.43 ± 0.02 
and i/p = 2.46 ± 0.03, which are also less accurate than 
1/ calculated with the critical Pc{G). Overall, Figs. 10 
and 11 illustrate the consistency of the RG approaches 
for the conduction distribution and for the level statis- 
tics, in the sense, that the best fixed point distribution 
of the level spacings corresponds to the fixed point of the 
conductance distribution. 



V. CONCLUSION 

Network models introduced in Ref. |5^ turned out to be 
a powerful tool to study the Anderson localization. With- 
out magnetic field, the propagation of electron waves 
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TABLE I: Part of fit results for v obtained from ai and cip 
for different system sizes A*', intervals around the transition, 
orders Ov and Oz of the fitting procedure. 



A [Zomin, 2:0max] Ov Oz V 



2 - 


512 


[9.93, 10.07] 


3 


2 


2.336 ±0.010 


2 - 


256 


[9.93, 10.07] 


2 


3 


2.412 ± 0.013 


4 - 


512 


[9.95, 10.05] 


3 


1 


2.325 ± 0.014 


2 - 


512 


[9.95, 10.05] 


2 


1 


2.402 ± 0.014 


2 - 


256 


[9.95, 10.05] 


2 


2 


2.360 ±0.016 


16 - 


512 


[9.95, 10.05] 


2 


3 


2.385 ± 0.018 


2 - 


128 


[9.93, 10.07] 


1 


3 


2.384 ±0.019 


4 - 


512 


[9.93, 10.07] 

ai 


2 


1 


2.471 ± 0.019 


2 - 


512 


[9.93, 10.07] 


2 


2 


2.383 ±0.010 


2 - 


512 


[9.93, 10.07] 


2 


3 


2.388 ±0.010 


2 - 


512 


[9.93, 10.07] 


3 


1 


2.346 ± 0.012 


8 - 


512 


[9.93, 10.07] 


2 


3 


2.376 ± 0.012 


2 - 


512 


[9.95, 10.05] 


2 


3 


2.368 ± 0.014 


2 - 


128 


[9.93, 10.07] 


2 


3 


2.377 ±0.016 


16 - 


512 


[9.95, 10.05] 


2 


1 


2.367 ±0.016 


2 - 


256 


[9.93, 10.07] 


3 


3 


2.372 ± 0.018 
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FIG. 11: Behavior of a.\ computed for initial distributions 
-PoaiiB different from the critical distributions, as shown in 
Fig. hol. Data are plotted for RG iterations n = 1, . . . , 9 cor- 
responding to effective system sizes A = 2" = 2,...,512. 
Curves for small n do not cross at the common point 20 = 0. 
Full lines indicate the functional dependence according to FSS 
using the minimization with Ov = 2 and Oz = 2. 
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FIG. 10: Comparison of the LSD Pc(s) and PGaufi(s) obtained 
from the corresponding conductance distributions shown in 
the inset. 

along each link of the network is allowed in both di- 
rections. In two dimensions the transmission coefficient 
of the network is zero fp& all parameters of the scatter- 
ing matrix at the nodes,E3 illustrating complete localiza- 
tion of electronic states. On the other hand, the two- 
channel network model with inter-channel mixing, that 
models spin-orbit inter^tion, exhibits a localization- 
delocalization transitionO thatis also in accord with the 
scaling theory of localization^ However, the version of 
the network model that has been most widely studied. 



is the chiral version, i.e., the CC model,c£l describing the 
electron motion in a disordered system in a strong mag- 
netic field limit. Within the CC model, the scattering 
matrix at the node is parametrized by a single number, 
e.g., the transmission coefficient t. On the qualitative 
level, the CC model yields a transparent explanation why 
delocalization occurs only at a single energy, for which 

— 1/2. On the quantitative level, in addition to the 
exponent, more delicate characteristics of the critical 
wave functions were exj-r acted from the numerical anal- 
ysis of the CC modelEj'Ell 

The fact that the RG approach, within which the cor- 
relations between different scales are neglected, describes 
the results of the large-scale simulations of the CC model 
so accurately, indicates that only a few spatial correla- 
tions within each scale are responsible for the critical 
characteristics of the quantum Hall transition. More pre- 
cisely, the structure of the eigenstates of a macroscopic 
sample at the transition can be predicted from the anal- 
ysis of a single RG unit consisting of only five nodes. 
Earlier we have dpnonstrated this fact for the conduc- 
tance distribution.L3 In the present paper this statement 
is reinforced by the study of the level statistics at the 
transition, which is a complimentary (to the conductance 
distribution) characteristics of the localization. 
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